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Abstract 

We construct a cosmological model which is physically reasonable, mathematically tractable, and extends the 
study of CDM models to the case where the equations of state (EoS) for matter and dark energy (DE) vary with 
time. It is based on the assumptions of (i) flatness, (ii) validity of general relativity, (iii) the presence of a DE 
component that varies between two asymptotic values, (iv) the matter of the universe smoothly evolves from an 
initial radiation stage - or a barotropic perfect fluid - to a phase where it behaves as cosmological dust at late 
times. The model approximates the CDM ones for small z but significantly differ from them for large z. We 
focus our attention on how the evolving EoS for matter and DE can modify the CDM paradigm. We discuss a 
number of physical scenarios. One of them includes, as a particular case, the so-called generalized Chaplygin gas 
models where DE evolves from non-relativistic dust. Another kind of models shows that the current accelerated 
expansion is compatible with a DE that behaves like pressureless dust at late times. We also find that a universe 
with variable DE can go from decelerated to accelerated expansion, and vice versa, several times. 
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1 Introduction 



One of the most challenging problems in cosmology today is to explain the observed late-time accelerated expansion 
of the universe [TJ-P3]. Since the gravity of both baryonic (ordinary) matter and radiation is attractive, the fact that 
the universe is presently accelerating, and may continue to do so, forces us to rethink and question some fundamental 
concepts about the universeQ 

The first question highlights our limited knowledge of the real nature of the content of the universe. In fact, an 
accelerated expansion requires the presence of a new form of matter (called dark energy), which could (i) produce 
gravitational repulsion, i.e., violate the strong energy condition; (ii) account for 70% of the total content of the 
universe; (iii) remain unclustered on all scales where gravitational clustering of ordinary matter is seen. (For a recent 
review see Ref. |18jV While we do not yet know exactly the physical mechanism responsible for this exotic behavior, 
we do have possible candidates for dark energy (DE). They include: a cosmological constant or a time dependent 
cosmological term jl9j-[20 ; an evolving scalar field known as quintessence (Q-mattcr) with a potential giving rise to 
negative pressure at the present epoch [2]- [25; dissipative fluids [23]; Chaplygin gas [25]- [U]; K-essence [2"5]-[3T]. 
and other more exotic models [32] . 

The second question is whether general relativity is applicable to describe the universe as a whole. Indeed, this 
and other puzzles of theoretical and experimental gravity have triggered a huge interest in alternative theories of 
gravity (See, e.g. [33J and references therein) where the cosmological acceleration is not provided by dark energy, but 
rather by a modification of the Friedman equation on large scales 34 , 35 . These include scalar-tensor theories of 
gravity [35J- [22], various versions of Kaluza-Klein theories, braneworld, STM and Brans-Dicke theory in 51? [33]-|46|. 

In view of these uncertainties, much of the research work being done in cosmology is based on the construction 
and study of specific cosmological models. The simplest one that predicts accelerated cosmic expansion and fits 
observational data reasonably well is the ACDM model [T]-[4] which is based on the assumptions of (i) flatness, (ii) 
validity of general relativity, (iii) the presence of a cosmological constant A and (iv) Cold Dark Matter (CDM). The 
main problem of this model is the huge difference between the observed value of the cosmological constant and the 
one predicted in quantum field theory. The other one, although not vital for the model, is that the assumption of 
CDM can not be applied to the entire evolution of the universe. 

In this work, we construct an alternative cosmological model where we keep the first two assumptions but relax 
the other two. Rather than choosing to investigate constraints on specific DE models, here we use a parameterization 
originally proposed by Hannestad and Mortsell [47] which can accommodate a number of DE models, including a 
cosmological constant. On the other hand, we employ a phenomcnological approach to describe the matter content 
of the universe as a mixture of different components, which are not required to expand adiabatically. The mixture 
smoothly evolves from an initial dense radiation stage - or a barotropic perfect fluid - to a phase where it behaves as 
cosmological dust at late times. 

This framework allows us to compare and contrast different physical settings and tackle a number of questions. 
Here we focus our attention on the evolution of the universe and its cosmic acceleration and on how the evolving 
EoS for matter and DE might modify the ACDM paradigm. We explore whether DE could play a significant role 
in the past evolution of our universe, and - conversely - whether the primordial EoS of matter can affect the details 
of accelerated expansion at late times. Also, whether dark energy can cross the line between quintessence and 
phantom regimes (the crossing of the cosmological constant boundary) [48j . Another important question is whether 
the accelerated expansion, once begun, continues forever: Is it possible for an ever-expanding DE-dominated universe 
to go through different cycles in which it changes from decelerated to accelerated expansion, and vice versa? 

This paper is organized as follows. In section 2 we give a brief introduction to Einstein's equations in a homo- 
geneous and isotropic background. In section 3 for the sake of generality we integrate the field equations without 
assuming any specific EoS. In this way we obtain general expressions for the Hubble, density and deceleration pa- 
rameters which in practice provide a simple recipe for the construction of cosmological models. In section 4 we 
introduce the EoS that generate our cosmological model and obtain explicit expressions for the relevant cosmological 

1 It should be mentioned that some authors keep a more skeptical point of view. They argue that the observational data, as it presently 
stands, can be explained without resorting to the existence of a negative-pressure fluid or a cosmological constant. The concept is that the 
departure of the observed universe from an Einstein-de Sitter model can be ascribed to other physical processes and/or to the influence 
of inhomogeneities. See e.g., | 15l . |16| . 1171 and references therein. 
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quantities. In section 5 we study the properties of a universe whose matter content is described by the EoS proposed 
in section 3, and the DE has a constant EoS. In section 6, in the context of CDM, we concentrate our study to the 
effects of a DE with variable EoS. In section 7 we present a summary of our work. 



2 Field equations 

A spatially homogeneous and isotropic universe is described by the FLRW line element, which in polar coordinates 
(r, 9, ip) is written atH 

^ ' r 2 (d9 2 + sin 2 9dp 2 



ds 2 = dt 2 ~a 2 {t) 



1 — kr 2 



(1) 



where a(t) is the scale factor with cosmic time t; k is the curvature signature which can, by a suitable scaling of r, 
be set equal to —1, or +1. 

The evolution of the scale factor and the stress energy in the universe are governed by the Einstein field equations: 

Rnv - -^g^R = T^v (2) 

The assumption of isotropy and homogeneity (J]} requires that the stress-energy tensor take on the perfect fluid 
forrro T{? = diagonal (p, — p, —p, —p), where p and p stand for the total energy density and total isotropic pressure 
of the cosmological "fluid" . This framework leads to two independent equations 

U , = f>«, (3) 



a 2 a 2 



1=1 



and 

2d 



N 

a a 



- = -P=-l^P {l) > W 



i=i 



where an over-dot indicates ordinary derivative with respect to t. The total energy density and pressure have been 
split up into constituents: and pW represent the energy density and pressure of the i-th component that fills the 
universe; the sums are over all, say TV, different species of matter present in the universe (baryonic, non-baryonic, 
radiation, cosmic neutrinos, dark energy, etc.) at a given epoch. 

These equations can be combined to obtain the continuity equation 

p + 3(p + p)-=0, (5) 
a 

which is equivalent to the covariant conservation equation T^f u = 0. 

Thus, there are two independent equations and (27V + 1) unknown quantities, namely, a(t), p«, pW, i = 1..N. 
To close the system we need to provide (27V — 1) additional equations. They could be TV equations of state (EoS) 
relating the pressures and densities. The remaining (TV — 1) equations are usually generated by the assumption that 
each component is expanding adiabatically, i.e., that there is no interaction between the cosmological constituents. 
What this means is that one imposes the energy conservation equation ([5]) on each constituent. 

The scale factor a(t) can be expressed as a Taylor series around the present time to: 

^ = 1 + H (t - t ) - \qoH 2 (t - to) 2 + \roH% (t - t Q ) 3 + (6) 
ao 2 b 

where H = a/a is the Hubble parameter, which is given directly by ([3]); q = —da/a 2 is the deceleration parameter 
that can be evaluated from (QJ. If the i-th component has an EoS p^ — ujWpW then 



2 In this work we use relativistic units: 8nG = c = 1. Also a subscript zero denotes the quantities given at the current epoch. 
3 To be rigorous, one should say that the matter has to be a fluid with bulk viscosity at most. 
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JV 



(0 



k 

2a 2 H 2 ' 



(7) 



where fiW = p(*)/3_ff 2 is the corresponding density parameter. Next, r = ot /aiJ is the statefmder parameter 
introduced by Sahni et al. [29] and Ulam et al. |50j . Taking time derivative in (j4j we obtain 



If, for the sake of generality, the interplay between the constituents is not neglected then each one evolves as 

JV 



1 + w 



p«=QW, with ^Q (j) =0, 



where measures the strength of the interaction of the i-th constituent with the rest of the components. Using 
this equation, the parameter r becomes 



JV 



'♦IE * 







dw^ ~) 




i + ™ w 


da J 



JV 



a 2 H 2 ' 



(8) 



In a similar way, one can express higher order terms in ^ as functions of w™, Q^' and their derivatives. 



3 General integration of the field equations 

The purpose of this section is to integrate the field equations with the least possible number of assumptions. With 
this in mind, it is convenient to split up the total energy density and pressure as 

P = 

p = p^+ P {M \ (9) 
where p^ DE ^ and p^ DE ^> represent the DE contribution and 

P w=x; P w > pw=^p« (io) 

i i 

are hybrids containing the contribution of relativistic particles, photons, the three neutrino species, as well as the 
contribution of non-relativistic particles (baryons, WIMPS, etc.). 

In term of these quantities the field equations ([3])- (|4j) now contain five unknowns. Here, we formulate no as- 
sumptions regarding the nature of the expansion of the constituents of the matter mixture p^ M ^; they may evolve 
non-adiabatically. However, to make the problem solvable, we neglect any matter-DE interaction and assume that 
the DE component is expanding adiabatically. As a result, p( DE ' and p DE , as well as the effective quantities p^' 1 ^ 
and p M , satisfy the continuity equation ([S]) separately. 

To close the system of equations we should provide two EoS. The simplest equation of state between density and 
pressure is the so called barotropic equation pj p = w, where w - in relativistic units - is a dimensionless constant. A 
direct generalization to this equation is to assume that w is not a constant but a function of the epoch. In this work 
we assume that both, matter and DE satisfy such type of EoS, viz., 

= w(a)p( M \ 

p( DE *> = W(a)p( DE \ (11) 
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To explain the accelerated expansion one has to accept that the DE component violates the strong energy condition^. 
Thus, in what follows we assume W(a) < —1/3. 

• With these assumptions the field equations ([3])-(|4|) become 

k 



DE) 



= 1 



a 2 H 2 



and 



?= 2 + 2 



2a 2 H 2 ' 



(12) 



(13) 



where tt^ = p^ /3H 2 and QS DE ) = p 1 - ^ /3H 2 . These can, formally, be regarded as two equations for ft( M ) and 
q(DE)_ Solving them we get 



2q—l — 3W k(l + 3W) 
3{w-W) ~ 3a 2 H 2 {w-Wy 
2q-l-3w k(l + 3w) 
3{w-W) + 3a 2 H 2 (w-W)' 



(14) 



We note that the denominator in these expressions is always positive because W < for DE, as well as for Chaplygin 
gas models. Thus, the fact that fi( M ) > and f^- - 6 ) > imposes an upper and lower limit on q, viz., 



1 + 3W 



k 



a 2 H 2 



<a< 



l + 3w 



k 



a 2 H 2 



In the epoch where < fiW < 1/2 the upper limit reduces to 

q A + \ [w + W] +2^H 2 



l + 3w {±) 



(15) 



(16) 



where «;M = W if k = 1, and w^- ) = w if k = — 1. For the ACDM modeH the first inequality gives — 1 < q < 1/2 
during the whole evolution and — \ < q < — 1/4 in the DE dominated era. 

• Given the EoS (fTTj) . the continuity equations for p( M > and p^ DE ' can be formally integrated to obtain the 
evolution of the energy densities as 



p (M) 
JDE) 



Ci - 3 f: 
—z-e J 



C 2 _ 3 rwM da 
— e J 



(17) 



where C\ and C2 are constants of integration. Let us use a* to denote the epoch at which p( M )(a*) = p( DB )(a*), 
and - for algebraic simplicity - introduce the dimensionless quantity 

x=— , (18) 



Using this notation, without loss of generality, we can write 



p (M) 

n( DE ) 



— e Ji 



C _ 3 P m* 



■ da 



(19) 



4 The strong energy condition for perfect fluids, in the comoving frame, requires p + p > 0, p + 3p > 0. 

5 In the ACDM model the universe is flat, filled with dust w = and the DE is attributed to the presence of a cosmological constant 
W = -1. 
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where C represents the common value shared by the densities at the transition poin10 x = 1 (a = a*) . Since W < 
and u> > 0, it follows that p( M ) > p 1 -- ^ for x < 1, and vice versa. 
In what follows we denote 



da. 



F(x) = e 3 /i ^ 
G(x) = e- 3 Ji — du . (20) 

• The more accepted interpretation of the observational data is that the current universe is very close to a spatially 
flat geometry (k = 0), which seems to be a natural consequence from inflation in the early universe. In accordance 
with this, in the rest of the paper we will consider a flat (k = 0) universe. 

The Friedmann equation ([3]) with k — becomes 

3H 2 = %r [F(x) + G{x)\ . (21) 
x 6 

Thus, the density parameters are 

q(M) = F ( X ) ry(DE) = G ( X ) (22) 

F(x) + G(x)' F{x) + G(x) { ' 

Since F(l) = G(l) = 1 it follows that (M) = Q( DE ) = 1/2 at x = 1, as expected. Now, taking the derivative we 
find 

d" (M) = 3[ff(x)-^)] „ Wn(i) g (23) 
dx x 

This shows that is a strictly decreasing function of x as long as [W^x) — w(x)] < 0, which is always the case 

for the mixture of DE (W < 0) and matter (w > 0). It goes from some 1/2 < f^ M '(0) < 1 to Q( M \x -> oo) > 0. 
What this means is that there is a one-to-one correspondence between fi^) and x. Therefore, if we know its value 
today then (|22|1 allows us to compute xq by solving the equation (as we mentioned before, the subscript zero denotes 
the quantities given at the current epoch) 



F(x ) 



n [M) 



n 



(DE) 



G(x ), (24) 



• The solution to this equation allows us to introduce the cosmological redshift z in term of which a = ao/(l + z). 
The relationship between x, defined in (|18j) . and z is 

z = — - 1. (25) 

a; 

In particular, the transition from a matter-dominated era to a DE-dominated one occurs at the redshift given by 

z* = x — 1. (26) 



Clearly, 



> 1, 




> 


for ti M) 


= 1, 




= 


for rt M) 


< 1, 




< 


for fi< M) 



< 1/2, 

-) =zo<! =1, z*=0 for^ M) = l/2, (27) 



6 For non-phantom (W > —1) and phantom (W < —1) models the equation p( M >(a f ) = p( DE \a,) has no more than one solution, 

because the DE energy density is strictly decreasing or increasing function of a, respectively [x dp ^ dx = —3(1 + W) p^ DE ^\. However, 
this is not necessarily so if the DE component has different regimes where it crosses the cosmological constant boundary (W = — 1). 
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Thus, xq > 1 indicates that the universe became DE dominated in the past; xq < 1 that DE will dominate in the 
future, and xq — 1 that fig^ = Qq DE ^ — 1/2, i.e., that the DE-domination begins right now. 
From (HI) and ([26]) we find 



dxo dz* Xq 



dn { DE) dn { DE) ~ 3n ( M) n { DE) [w(x ) - w(x )] 



(28) 



which is always positive. This shows that the closer is fig 13 ^ to 1 , the higher the redshift z* - i.e., the further back 
in time is the transition to the DE-dominated era. 

Now we can evaluate the constant C. Using (|2"TT) and (|2"4")1 we find 

■m 2 x 3 n (M) iH 2 x 3 n (DE) 

F(x ) G(x ) ■ (29) 



With this quantity at hand, the Hubble parameter (|2"Tj) can be written as 



'•( 



2^0 \ 3/2 



H = H L . 

x / 



(M) F(x) n (DE)G{x) 
U ° F(x ) +S2 ° G(x ) 



1/2 

(30) 



• The deceleration parameter can be obtained from this expression by using 



x dH „ . . . 

"=-H^-^ (31) 



or directly from (|13p with k = 0. Either way we find 



[l + 3w(x)}F(x) + [l + 3W(x)}G(x) 
q{x) = 2[F(x) + G(x)} ' (32) 



In a similar manner, the deceleration parameter today is obtained immediately from (|13p with k = 0, as 

q = ~ {l + 3 [n^wixo) + ni DE) W(x )] } , (33) 

where xq is the solution to for a given £Iq M ^ . 

The function i^(x) decreases, and G(x) increases, with the increase of x. Besides W < —1/3. Therefore, at some 
point q changes its sign from positive to negative. Let us use x to denote the root(s) of the equation 

q(x) = [1 + 3w{x)} F{x) + [1 + 3W(x)] G(x) = 0. (34) 

It should be noted that it may have several real roots because q is not necessarily a monotonically decreasing function 
of a;. In fact, from ([32]) we get 

d q 3r>( M ) dw 9n^n^ . TTn3 30( DE ) dw 

t = — t. ; ^ \w—W\ H — (35) 

dx 2 dx 2x L J 2 dx 

In this expression the first term is non positive, (dw/dx) < 0, since the EoS of ordinary matter should become softer - 
or at most remain constant - during the expansion of the universe. The second one is negative. The third term could, 
in principle, be positive, zero or negative depending on the specific model. Therefore, an EoS for DE having a large 
positive slope may lead to several real solutions in (|34[) and the universe would pass from deceleration to acceleration, 
and vice versa, several times during its evolution. In general, the only sure thing about the global behavior of q is 
provided by the inequalities ()15[) and ()16[) . In summary, if (dW/dx) < then q is a strictly decreasing function of x 
and one can affirm that (fMf has only one solution. 
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Finally, we note that the acceleration of the universe at the end of the matter dominated phase, is given bj0 

q* = g(o.) = q(x = 1) = ~ + Hi) + W(l)] , (36) 

where we have used that F(l) = G(l) = 1. Clearly, if < then the universe comes into an accelerated expansion, 
at some X < 1, during the period when matter still dominates over the dark energy component. In contrast, if > 
then the accelerated expansion begins, at some x > 1, only after the DE component dominates over matter. If g» = 
then the onset of accelerated expansion coincides with the end of the matter-dominated era. 



4 Cosmological Model with variable EoS for matter and DE 

The formulae developed in the preceding sections give the cosmological parameters (iJ, tt, q, r) explicitly in terms of 
w, W, F and G. Also (|24l) -(|25 |) allow to express them as functions of the redshift. What we need now, to construct a 
cosmological model, is to have the appropriate EoS. In this section we do three things. First, we introduce a variable 
EoS for the matter mixture (jlOl) . and discuss its features. Second, we study a simple parameterization for DE where 
W varies between two asymptotic values. Third, we present the explicit expressions for the Hubble, density and 
deceleration parameters. 



4.1 Effective EoS for the matter mixture 

In the approximation where relativistic matter is modeled by radiation and non-relativistic matter by dust, to 
describe the evolution of the universe from an early radiation-dominated phase to the recent DE-dominated phase 
the function w must satisfy the conditions w — > 1/3 for and w — > for x 3> 1. Certainly, there are an infinite 

number of smooth functions that satisfy these conditions. 

To motivate our parameterization of w let us go back to (jlOp . The two main components of the cosmic mixture 
are provided by relativistic and non-relativistic particles. Therefore, 

J2 ■ pW p (NR) + p (R) 

w = ~ p (nr) + p {r) ' 

where p^ R \ p^ and p( NR \ p( NR } denote the energy density and pressure contributed by relativistic and non- 
relativistic particles, respectively. In the approximation under consideration = 3p( R ^ and p( NR ' = 0. Thus, 

1 „ s _ P [NR) 

W ~3[/(.*) + l]' m pW 

The relative contributions of the different components depend on time. During adiabatic expansion p( NR > decreases 
as x -3 , while p^ does it as x~ 4 . Thus, f(x) ~ x - at least - if one ignores interactions where an individual kind 
of particle becomes non-relativistic, gets bound, or annihilates. In this paper we assume f(x) — Kx a , where K 
is a constant of proportionality and a is a positive constant parameter. Thus, we adopt the following very simple 
approximation 

W = 3(^ + 1)- (37) 

The parameter a determines the rapidity of the transition from radiation to dust as well as the duration of the 
radiation-dominated epoch. In fact, differentiating with respect to x we obtain 

1 dw _ aKx^- 1 ^ 
w dx ~ Kx a + 1 



7 It is interesting to note that q„ depends only on the specific values of the EoS at x = 1 and not on F and G, which contain the 
information on the overall effect of the EoS on the evolution of universe. 
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Figure 1: The effective EoS (|37|) for the matter mixture (JTOj) . The Figure shows that for large values of a the 
transition from radiation to dust is very abrupt, and vice versa. The EoS at x = 1 is given by w, = 1/3 (K + 1). 
For the sake of illustration, the graphic was constructed with K = 10. 



What this means is that w has a plateau, namely w ~ 1/3, in the early universe for any a > 1. For a < 1 there is 
a rapid variation of w at x <C 1 - there is no radiation-dominated plateau - and the transition to dust occurs very 
slowly. Besides, the larger the choice of a > 1 the faster the transition to the dust era w ~ 0, which goes on with 
the increase of x. See Fig. [T] 

The coefficient K determines the EoS at x = 1, viz., 

w * = w{x = 1) = HihTy 

To have lo.siO we assume that K is a large, but finite, number. 

The mathematical simplicity of (|37|) allows us to express all the relevant physical quantities in terms of elementary 
functions. In fact, substituting (|37l) into (fT9|) we find 



Consequently, 



p ( M) = Ci(^°' + l)' / ° i „ #0| Cl = (K + l)-"°C. (39) 



p(M) = Cl( ^ 3+ ;) < -° v ° . («) 



• As expected, for a = 1 the constituents of the matter mixture recover their individual identities and 
separates into radiation and dust in adiabatic expansion, viz., 



For any a/1, the EoS giving the relation between p( M ) and /t/ M ) is given in parametric form by (|39|) and ((40 
The explicit expression p( M ) = p( M ) (p( M >) is quite cumbersome. However, the asymptotic behavior is as follows: 
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• For 



• For x > 1, 



JM) _ ^1 
—A ' 



„( M ) ~ i.„( M ) 



C x #(l-a)/a 

3a;( 3 + Q ) : 



For a > 0, the matter pressure decreases with the expansion of the Universe faster than the density. Note that p^' 1 ^ 
decreases as in cosmological dust models. However, formally, exact dust models correspond to the limit a —¥ oo. 
Bellow, in Section 5.1, we will discuss the properties of the models having a>l. 
• For i»l, near the transition density p* the EoS can be written atH 



1M) 



(3K + 4) [p^] 



(a+l) 



3(#+l){2(# + 2) [pW]' 



A' 



"} 



This EoS closely mimics the behavior of the one discussed by Israelit and Rosen 
density vary continuously through different epochs in FRW universe models. 



(42) 

where the pressure and 



4.2 A variable equation of state for DE 

To describe the evolution of DE we adopt the EoS proposed by Hannestad and Mortsell [47] , which in our notation 
become^ 



W 



U!X P + 7 

+ 1 ' 



(43) 



where /?, to and 7 are constants parameters. If w = 7 then W = u>. In any other case ui and 7 describe the asymptotic 
behavior of W. Namely, for f3 > 

W-J W !° rlM ' (44) 
[7 for a: « 1. v y 

For /3 < 0, the constants w and 7 interchange their role, i.e., W — > 7 for x 3> 1 and — > u) for a; -C 1. In what 
follows without loss of generality we assume (3 > 0. 
Besides, 

dW _ fix?- 1 (uj - 7) 
~dx~ ~~ ( x f> + l) 2 

shows that for /3 > 1, VF undergoes a rapid transition from an early epoch where W ps 7 to a large epoch where 
VF « w. In addition, W is an increasing function of x if > 7, and vice versa. 

Next, to account for the repulsive nature of DE we should demand the violation of the strong energy condition of 
classical cosmology. It is easy to verify that p^ DE ^ + 3p( DE ^ = p 1 - ^^! + 3W) < in the whole range of x provided 



1 

W< -3' 



1 

7< -3- 



(45) 



8 From (39} it follows that y 4 - (Ci/pO 1 ') " Ky - (Ci/p (M) ) ° = with y = x a . At the transition poiny x = 1, y = 1, ^Ci/p!; 1 ' 

1/ (if + 1). In the neighborhood of y = 1 we can set (Ci/p( M ') =1/ (X + and y = where £ is a dimensionless small param- 

eter. Substituting into the equation for y we get j/i = (if + l) 2 / (3if + 4), to first order in £. Now £ = [(Ci/p' M ') — 1/ (if + 1)1 and 



using that (Ci/p 



(Ci/p( M )) a = (if + (pi M) /P (M) ) Q we obtain y = x a = (3if + 4)" 1 (2if + 3) + (if + 1) (pi M, /P (M) ) 



Finally, 



replacing this into l|37j l and using JTTJ we get 02}. 

9 In our notation the four constants wo, wi, q and a B in Hannestad- Mortsell parameterization I47| can be expressed as q = /3, wq 

u>1 = 7, a? = (lUi/lOo) a* • 
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Lower limits on these quantities, namely u> > — 1 and 7 > — 1 are obtained if one assumes that DE satisfies the 
dominant energy condition p( DE ^ > |y P ' E l : 
Substituting (g3J into CEH]), we obtain^] 



P {DE) - 2 ( ^(i +7 ) . 0*0, C 2 = 2 3 (^C. 



Thus, 



P 



C 2 (ra^ + 7) 



2,3(1+7) ( x /9 + ^[3(^-7)+^]/^ ■ 

We immediately notice that for j3 — 3 (7 — uS) the above equations reduce to 

1 1 



(46) 
(47) 



P {DE) = C 2 
(DE) = c 



P 



x 3(l+u) x 3(l+7) 
U! 7 



(48) 



x 3(l+w) a;3(l+7). 

What this means is that for /3 — 3 (7 — u>) the DE can be interpreted as the superposition of two fluids with EoS 

pi=upx, p 2 =7P2, (49) 

each of which satisfies the continuity equation. 

The asymptotic behavior of (|46|) - (l47p is as follows: 

• For i<1, there are two different physical situations. If 7 7^ we get 



P 



(DE) 



C 2 



2,3(1+7) ' 



In this limit p( DE ^ / >M w x 1 "^. Therefore p< DS ) < p^ for x < 1 as long as 7 < 1/3. 
If 7 = we find 



(DE) ^ C2 



Thus, for j3 > 3 the DE component - in this limit - behaves exactly as in cosmological dust models. 
• For i>lwe obtain 



JDE) 



Co 



a;3(i+u) • 



Since u < 0, it follows that > p( M ) and |p(-° B )| > p( M ) for a; > 1. Also, in this limit when u = -1 the DE 

component behaves as a cosmological constant p( DB ) = —p^ DE ^ = C2. 

• For the transition from VF ~ 7 to VF ps u> can be approximated by an EoS similar to (1421) . 



Generalized Chaplygin gas: The above analysis indicates that when lj = — 1, /3 > 3 and 7 = the DE component 
evolves from 11011- relativistic matter (dust) to a cosmological constant, similar to the so-called Chaplygin gat0. In 
fact, it is not difficult to see that for this choice of parameters P5|) - (jT71) yield 



SDE) 



= -a 



(50) 



10 The choice = correspond to a constant equation of state W = (uj + 7) /2. 

11 Chaplygin gas satisfies the EoS equation of state p = —A/p" where p is the pressure, p is the density, with a = 1 and A a positive 



constant. In generalized Chaplygin gas models a is a parameter which can take on values < cr < 1 
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where cr = (f3 — 3)/3 and 



Thus, the EoS 



P 



W = C 2 



1 



w = 



X? 



i/(i+^ 



(51) 
(52) 



with 3 < (3 < 6 describes a generalized Chapligyn gas. 



4.3 The Hubble, density and deceleration parameters 

We now proceed to apply our general formulae to the EoS (|5T]) and (113")) . For practical reasons, it is useful to write 
(157)1 as 

w = - r -, (53) 

[Kx a + 1) ' v ; 

with n = 1/3. This will allow us to follow in detail, step by step, the possible effects of the primordial radiation 
on the observed accelerated expansion today. Besides, it allows us to extend our results to include other types of 
matter, e.g., stiff (incompressible) matter (n — 1). Also for n = we should recover the usual CDM picture. 
Thus, for the sake of generality in our discussion we let 



< n < 1, 

in accordance with the strong and dominant energy conditions. 
From (gOJ) we get 



(54) 



F(x) = (K + iy 3n/c 



K 



1 



-\3n/a 2 3 ( aJ ~f)/' 3 



-3(w- 7 )//9 



(55) 



• The Hubble parameter is obtained by substituting these into (|30[) . namely, 



H(x) = H 



^ M) (^) 3 /w+^r ) (^) 3(w+1) ^) 



1/2 



(56) 



where we have introduced the functions f(x) and g(x) defined by 



/(*) 



K 



^ -i 3n/c 



/(*) 



X 

1 + 



3(w-7)//3 



(57) 



and xq is the solution to (j2"4"]l , which in the case under consideration becomes 



x 



3n/a 



3(w- 7 )//3 



1 + — a 



_ 2 3(w-7)//9 f K + ^3n/a 



n 



(M) 



L n 



(DE) 



(58) 



The energy densities, from ([T9| . ff20|) . ([55]) . ([57]). can be written as 
For n — and 7 = w, these expressions reduce to the usual ones in the CDM model, as expected. 



(59) 
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Figure 2: The figure illustrates the situation where the deceleration parameter changes sign several times during the 
expansion of the universe. We have selected ui = —1/2, 7 = —3/2 and /3 = 20 to generate a large dW/dx > 0. Also, 
we have chosen a = 100 to assure that the models with primordial radiation (n = 1/3) and stiff matter (n = 1) are 
indistinguishable from the dust model at late times. Also, for this choice the number of zeroes of q(x) as well as 
their locations are essentially independent of the specific value of K. The graphic was constructed with K = 1. 



The deceleration parameter is easily obtained by a direct substitution of (|43l) . (|53"f and (|55j) into (|3"2"j) . We get 

3w ( K , J_\3n/a l~i^*+i+3„"| ^S^-tV* 3 (if + 1 ) 3n / a (\ + 1 ) ~^-~tW [ (3n+l)s g + (37+1 
\ ' x a I Kx a + 1 ' \ ' / V ' ~xfi I xf + 1 

= r— r-w^ — : r^zzzvm • ( 60 ) 



x ^ (^ + _L) 3 «A* 


Ksc a +l+3n 
Kx a + l 


+ 2 3 ("-^)/' 3 (X + 1) 3 " /Q (1 + ^)- 3( "~ 7)//3 


"(3w+l)a:' 3 + (37+l)" 


2 


>. + _L) 3 ™/« + 2 3( W -7)//3 + l) 3 "/« (1 + J^) 





We observe that q(x = 0) = (1 + 3n) /2 and q(x —> 00) = (1 + 3w) /2. Since w < —1/3, there is at least one value 
of x, say x — x, for which q vanisheJ^l. In accordance with our discussion in (j35j), q is not necessarily a decreasing 
function of x if dlF/cfcc > in which case q{x) = may have more than one real root. In Figure [2] we illustrate the 
situation where q(x) = has three different roots for different EoS: dust (n = 0), primordial radiation (n — 1/3) and 
stiff matter (n = 1). 

In general, only the upper and lower bounds on q are well-defined and arc given by (1151) . In the case under 
consideration these are 

<g<l, £ = min{7,Lj}. (61) 



Besides, in agreement with (|36[) we find that at the transition point x — 1 (a = a*) 

1 3 (uj + 7) 3n 
M(;C = 1)= 2 + ^^ + 4WTl)' 
regardless of the specific choice of a and f3. 

• The magnitude of the deceleration parameter today, qo, is readily given by (]33[). viz., 



(62) 



1 

go =2 



3nf2 



(M) 30 



(DE) 



(wx$ + 7 



(63) 



2 In Section 5.3 we will consider models with ui = 0. 
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The same expression can be obtained from (I60p . after using (|58l) and some algebraic manipulations. 



4.3.1 Cosmological parameters in terms of the redshift 



For practical/observational reasons it is convenient to express the cosmological parameters as functions of the redshift 
z introduced in ([25]) . 

The Hubble parameter (|56|) becomes 



H 2 = H$ 



with 



M 

9(z) 



Kx% + (1 + z) c 



Kx% + 1 
x$ + (1 + zf 



3n/a 



4 + \ 



-3(w- 7 )//3 



The energy densities can be written as 

P 

Finally, the deceleration parameter as a function of z is given by 



P 

SDE) 



3(l + zf4 M) H*f(z), 

3(i + ^nrffo 2 s(4 



(64) 



(65) 



(66) 



(1 + 3w) + (1 + zf (1 + 37) 



a£ + (l + *f 



(67) 

These equations show that at low redshift /(z) w 1 and g(z) ps 1 and we recover the usual CDM model, i.e., the 
dynamics of the EoS for DE and matter only becomes important at hight redshift. 



5 Universe of matter with variable EoS and DE with W = constant 

We now start the study of the properties of our cosmological model. In this section we concentrate our attention on 
the question of how a dynamical EoS for matter can alter the details of the accelerated expansion as given by the 
CDM (w — 0) models. With this in mind, here we confine the discussion to the case where the DE component has 
a constant EoS (a variable W will be discussed in the next section). 
Thus, in this section 

W = oj = constant. (68) 
For ui = — 1 we have a standard cosmological constant. Setting 7 = ui, from ([5H|1. (foT))) and (f6"3"| we get 



3n / a 



3n/a ^0 



(M) 



(DE) ' 







q(x) = x 3 



K 



3n/a ^ Kx a + l + 3n 
Kx a + 1 



+ (K + lf n/a (1 + 3w) = 0, 



(69) 
(70) 



and 
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1 r 



i + 3ujn { DE 



(M) 



2{Kxf ) + 1)' 



(71) 



respectively. Thus, to obtain a numerical value for qo as well as the observational quantities z* = xq — 1 and 
z = xq/x— 1, we need to solve and ([TO)) . 



CDM: For n — 0, they have a close algebraic solution which is the CDM model, viz. 



(CDM) 



''() 



(M) 



-i l/3w 



o 



(CDM) = [_(l + 3w )]V&- 
1 

2 



(CDM) 
% 



1 + 3wfi^ DB) 



(72) 



Variable EoS: For n ^ 0, from (|M)) - ([7D|) we find expressions for xo and a; which are very similar to those given 
by (|72j) . namely 



x = 



»'o 



n I a(ui — n) 



n/ct(u)—n) 



(M) 



n 



(DE) 



l/3(u-n) 



1 +3w 
l + 3u> 



l/3(w-n) 



(73) 



where u>o, u>* and u> represent the EoS (|53|) evaluated at Xo, a; = 1 and a; = x, respectively. They reduce to their 
CDM counterpart for n = 0, as expected. Although the shape of these expressions is intriguing, to find xq and x we 
need to provide a set of values (S1,q M \ u), n, K, a) and then solve numerically. 

However, a direct inspection of (pi?)) - (1701) suggests that we examine the following cases: (i) a = 3n; (ii) < a <C 3n 
and (hi) a 3> 3n. 

• When a = 3n, the matter in the universe is a superposition of two non-interacting fluids; one of them is dust, 
and the other one is a fluid with an EoS p = np. Specifically, 



,(M) _ 



KCi d 



+ 



When < a <C 3n we have 



X 3 a;3(ri+l) ' 

Cl 



P 



(M) _ 



nCi 

x 3(n+l) ' 



(74) 



x 3(l+n) ' 



K + l 



x 



n 



(M) 



(DE) 
L 12 



l/3(w-n) 



-(l + 3w) 



-i l/3(w-fi) 



1 + 3n 



(75) 



(76) 



Since if is a large number (X ~ p( NR )/p( RS > ~ 3 x 10 4 today), the above solution is practically indistinguishable 
from the CDM model (f72|) . 

• When a 3> 3n, we find that the physics is independent of the particular choice of K and a. In the next 
subsection we discuss the details. 
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Figure 3: The deceleration parameter q versus x for models with cosmological constant. The figure compares the 
evolution of q in the ACDM model (n = 0) to the one in models with n = 1/3 and n = 1 calculated for a = 100 and 
K = 1. 



5.1 Asymptotic model: a > 3n 

A simple numerical analysis of (|69 |) - ([70| . with (Qq M \ uj, n, K) fixed, shows that xq and x tend to specific finite 
limits when a — » oo. A further analysis demonstrates that these limits are insensitive 13 ! to the specific choice of K. 

Therefore, in this limit only the parameters uj, n and Qq M ^ have physical relevance. It corresponds to the case 
where the transition to dust occurs abruptly near x = 1 (See Figures [TJ [3j |4] ) . 

As shown in (|27j) . the character of the solutions to (|69j) strongly depends on whether the universe is dominated 
by matter or DE. Accordingly, for a>l we find 



(CDM) 



for Q 



(M) 



o < 1/2, 



for ft 



(M) 





> 1/2. 



(77) 



For fig*^ = 1/2, xo — 1 regardless of the choice of parameters. To simplify the discussion in Table 1 we calculate xq 
for various < 1/2 and u. 



Table 1. xq for various ft M ', w, any n and a 3> 1 




w = -0.4 


u = -0.5 


uj = -2/3 


w = -l 


uj = -1.3 


uj = -1.5 


uj = -10 


0.1 


6.240 


4.327 


3.000 


2.080 


1.756 


1.629 


1.076 


0.2 


3.175 


2.520 


2.000 


1.587 


1.427 


1.361 


1.047 


0.3 


2.026 


1.759 


1.527 


1.326 


1.243 


1.207 


1.028 


0.4 


1.402 


1.310 


1.225 


1.148 


1.110 


1.094 


1.014 



• By virtue of ([35]). the EoS (|68|) guarantees that q{x) is a strictly decreasing function of x. Consequently, ([70]) has 
only one real root for each set (uj, n, K, a). They heavily depend on whether the deceleration-acceleration transition 



13 In fact, one can show that J- feft ~ IM, 

1(1 dK X aK 



so that -S~i#^0asa^oo. 

xn aK x aK 
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Figure 4: The statefinder parameter r, given by (jSJ), versus x for models with cosmological constant. In the ACDM 
models r — 1. The sharp peak of r is evidence of the rapid change experienced by the deceleration parameter. Once 
again we have taken a = 100 to ensure that the models behave like dust at late times. Moreover K = 1. 



occurs before or after the end of the matter-dominated era. From (l36l) we find 



4 V K + l 

What this means is that x > 1 for to > — [2/3 + n/ (K + 1)], and vice versa. 
Thus, for a 3> 1 we obtain the following solutions: 

(i) u > -2/3: 

.f = [-(l + 3w)] 1/3 "\ 
where the condition u > —2/3 ensures that x > 1. 

(ii) - [2/3 + n/(if + 1)] < w < -2/3: 

x = 1+. 

(iii) w = - [2/3 + n/(JC+l)]: 

x = 1, 

for any choice of parameters. 

(iv) - (2 + 3n) /3 < co < - [2/3 + n/(K + 1)]: 

x = 1". 



1 / 3n , 

? „ = - 2 + 3W+-— - . (78) 



17 



(v) cj < - (2 + 3n) /3: 



(l + 3o;)' 
(l + 3n). 



l/3(u-n) 



where the condition w < — (2 + 3n) /3 ensures that x < 1. 

For n = we recover the expressions for CDM. In Table 2 we present the approximate values of x for different 
choices of cj, n and a 1. We remark that those for w > —2/3 are independent of n and coincide with the CDM 
values. 



Table 2. a; for various cj, n and q»1 


w 


n = 


n = 0.1 


n = 0.2 


n = 1/3 


n = 2/3 


n = 1 


-0.4 


3.824 


3.824 


3.824 


3.824 


3.824 


3.824 


-0.5 


1.587 


1.587 


1.587 


1.587 


1.587 


1.587 


-2/3 


1 


1 


1 


1 


1 


1 


-5/6 


0.850 


0.950 


1 


1 


1 


1 


-1 


0.794 


0.878 


0.940 


1 


1 


1 


-1.3 


0.761 


0.826 


0.876 


0.927 


1 


1 


-1.5 


0.757 


0.813 


0.858 


0.903 


0.976 


1 


-10 


0.894 


0.903 


0.910 


0.917 


0.932 


0.942 



If cj < —2/3, then our model predicts some changes to the CDM paradigm. First we note that in CDM models 
(first column in Table 2) x varies with cj. But, this is not necessarily so when n ^ 0. Solutions (ii)-(iv) show that 
x « 1, for every n, in the whole range — (2 + 3n) /3 < cj < —2/3. 

However, the redshift associated with x does depend on cj through xq. For example, if we take n = 1/3 and 
Qq = 0.3 then Table 1 gives xq = (1.527, 1.326) for u = —2/3 and cu — —1, respectively. Correspondingly, we get 
z = (0.527, 0.326), although x = 1 in both cases. 

When cj < — (2 + 3n) /3, for Qq < 1/2 the redshift for the deceleration-acceleration transition is given b}Q 





l/3w 


(1 + 3cj)" 


l/3(n— lj 


o( £,£ ) 




(l + 3n). 





For the case where DE is a cosmological constant (w = —1) and f2 = 0.3 we get 



z = (0.671, 0.511, 0.411, 0.326) , for n = (0, 0.1, 0.2, 1/3) . 

We note that the redshift at q = predicted by the CDM models is a little more than twice the one obtained for 
n = 1/3. As a second example we consider phantom DE with cj = —1.5 and, once again, take = 0.3. We find 

z= (0.595, 0.484, 0.407, 0.337, 0.237, 0.207), for n= (0, 0.1, 0.2, 1/3, 2/3, 1). 

Thus, for non-phantom [cj > —1) and phantom matter [cj < —1), we can conclude that the stiffer the primordial 
EoS, the later - i.e. closer to our era - the transition to accelerated expansion. 

To finish the discussion, we would like to point out that the asymptotic model discussed here and the CDM model 
- given by (|72|) - do not challenge each other. Instead, they are complementary. Indeed, they represent two different 
limits: in the CDM model dust is assumed during the whole evolution, while here there is a sharp transition from a 
primordial stage, to cosmological dust. These two asymptotic models should serve to restrict or restrain more 

realistic ones. 

14 For n ( Q M) > 1/2 we should use the second solution in ((23 • See also (188 D . 
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5.2 Features of the model that are independent of a and K 

In any physical model is important to identify the features that are independent of the particular choice of the 
parameters of the theory. Therefore, we now look for relationships between the observational parameters that are 
independent of the particular choice of a and K. With this in mind we recast (1691) into the form 



X 



3u 



>4 



-i 3n 



K 



?>UJ 



>4 



3/ !, 



a 



(M) 







(DE) 



where X — Kx a . Since, a > 0, n > and u) < 0, from this equation it follows that X > K for Oq ^ < n Q DE ^; 

X < K for Q M) > fl { DE \ and X = K for O M) = fl { DE) . This is consistent with the fact that x > 1 (x < 1) 
0^ < 1/2 (Off > 1/2) discussed in {27}. 
Now, from (|7T|) we obtain (n 7^ 0) 



3 (u - rc) O M) + 2g - 1 - 3w 



3wO, 



(DE) 



2?0 + l 



(80) 



Using this expression we find: 



Models with o[, M) < 1/2: Setting ifx > if, from {80]) we get 



1 + 3wO 



< % < 



1 



l + 3wO 



(DE) 



-nO, 



(M) 



(81) 



As in the CDM models, the current accelerated expansion requires 

1 



uj <-- 



30 



(DE)' 



(82) 



Equivalently, in terms of uj the above inequality can be written as 



2gp-l 



raO 



(M) 



20, 



< cj < 



2«lb-l 



30 



(DE) 



(83) 



For any given Oq*^ < 1/2, this expression yields a lower and/or an upper bound on qo for various values of u, 



-1 + 



30 



(M) 



for ui > — 1, 



3O M > 



(M) 



9o < 



< <7o < -1 



3 (n + 2) ( 

4 ' 
3 (n + 2) O M) 



for a; < — 1, 
for ui = — 1. 



(84) 



As an illustration, let us take Oq ^ = 0.3. Then, accelerated expansion requires u> < —0.476; if U) > —1 then 
the deceleration parameter has a lower bound q^ 71 ^ — —0.550, regardless of n; if uj < —1 it has an upper bound 

that depends on n, namely, q^ nax%> — (—0.325, —0.447, —0.550) for n = (1, 1/3, 0), respectively; if lj = — 1 then 
qo is bounded from above and from bellow, viz., —0.550 < qo < (—0.325, —0.447, —0.550) for n = (1, 1/3, 0), 
respectively. 
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Models with U 



(M) 




1/2: In these models the condition xq — 1 yields 



1 + - (n + 2oS) 
4 



Accelerated expansion at the current epoch requires 



to < 



4 + 3n 
6 ' 



(85) 



(86) 



Models with O,^ 1 ^ > 1/2: In these models the lower and upper bounds are generated by the condition < x a < 1. 
Using (JHOJ) we finq^l 



1 + 3iun { DE) 

)(M) 



— nfl 



(M) 



<qo< 2 



1 + 3wft 



(DE) 



—nfl 



(M) 



(87) 



This expression, together with S7 > 1/2, can be used to obtain more specific bounds on qo, similar to those in 

©. 

Although observations indicate that the universe today is DE-dominated, it is of theoretical interest to note that 
accelerated expansion may still occur in a matter-dominated phase if |w| is large enough. Specifically, 



UJ < - 



in Q 



(M) 



m[ DE) 



(88) 



As an illustration, for Qq^ = 0.6 the transition from decelerated to accelerated expansion occurs in a recent past 
if u < - (5/6 + 3n/4) = - (0.833, 1.083, 1.583), for n = (0, 1/3, 1), respectively. If u > - (5/6 + 3n/4), due to the 
continuous expansion of the universe, the transition occurs in the future, i.e., at some —1 < z < 0. In the above 
expressions setting n = 0, we recover the CDM model. 



6 Dark energy with variable EoS and CDM (w = 0) 

In this section we consider the case where the DE component evolves according to the EoS (|43[) . Our aim is to probe 
the extent to which a dynamical EoS for DE can affect and/or modify the straightforward description of accelerated 
expansion as given by the W = uj = constant models. As a framework for our discussion we consider a flat universe 
whose matter content is approximated by non-relativistic dust. 

Thus, in this section we set n = 0, and keep uj 7^ 7. With this simplification (|58|) reduces to 




3(w-7)//3 



2 3(w-7)//3 



(M) 



(DE) 



n 



From (1601) . with n = 0, we get the equation for x, viz., 

-3(«-7)/fl r(i + 3w )^ + (i + 37 ) 



-3w , 23(w-7)//3 



1 



1 

w 



0. 



(89) 



(90) 



For 7 = uj we recover the CDM model ((72)) . 

As we mentioned in Section (3.2), the EoS (T43l - although in a slightly different notation - has been studied by 
Hannestad and Mortsell in Ref. [47] . From a joint analysis of data from the cosmic microwave background, large 



15 x > implies 1 + SuiQ^ 3 ^ < 2q < 1 + 3wQ<, DB) + 3nf^ M) 
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scale structure and type-la supernovae, these authors evaluated the constraints on f|43[) and concluded that the best 
fit model corresponds, in our notation, to u) = —1.8, 7 = —0.4, /3 = 3.41 and fig = 0.38. For these values from 



-(|90p we get Xn = 1.148 and x = 0.789. Accordingly, the universe becomes DE-dominated at z + = 0.148 and 
starts accelerating expansion at z = 0.452. From ((53"]) we find q$ = —0.673 and Wo = —1.262. 

However, the observational facts are not yet conclusive as to settle the question of the dynamics in the dark 
energy EoS. Indeed there are other choices of parameters that lead to results consistent with observational constrains. 
Therefore, here we concentrate our attention on studying the different possibilities and models for DE offered and/or 
suggested by (|43|) . 



A short inspection of (l89|) - (l90l) suggests we start with the examination of three different cases: (i) /3 = 3 (7 — w); 
(ii) < (3 -c 3 (7 — uj) and (iii) (3 ^> 3 (7 — uj). The first case corresponds to the superposition of fluids mentioned in 
(|48p. The second one yields the CDM solution (|72p with uj replaced by (uj + 7) /2. The third case characterizes the 
limiting case where W rapidly evolves from 7 to w. 

6.1 Quick transition from W = 7 to W = oj 

We now proceed to study the properties of the limiting model. For /3 > 1 the solution to is given by 

l/3u> 



x 





' n (M) ' 




(DE) 
_ S2 




" n (M) ' 




n (DE) 
_ s2 



1/37 



for O M) < 1/2, 
for Q M) > 1/2, 



(91) 



where we have assumed 7 / and w^0. The solutions with 7 = and uj = will be discussed bellow. 

Regarding (|9"0)) , for /3 ^> 1 we find several solutions depending on whether uj and 7 are greater or less than —2/3. 
These are0 : 

(i) For uj > —2/3, 7 > —2/3 the solution is x = [— (1 + Sw)] 1 ^ 3 ^. The requirement uj > —2/3 asserts that x > 1. 

(ii) For w > —2/3, 7 < —2/3 there are three distinct roots: x\ — [— (1 + 37)] 1 / 37 ; one of the following x-z — 
(1 + , 1, 1 _ ), depending on whether g* < 0, = or > 0, respectively; and £3 = [— (1 + Sw)] 1 ^ 3 ". The conditions 
on uj and 7 establish that xi < 1, X3 > 1. Also, here 

1 3 

q* = 9(1) = 2 + g ( w + 7) • 

Thus, for a given oj the statement > denotes 7 > — (uj + 4/3), and vice versa. 

(iii) For oj < —2/3, 7 > —2/3 we find either x — 1 + or x — 1" depending on whether > or <?* < 0, respectively. 

(iv) For uj < —2/3, 7 < —2/3 the solution is x = [— (1 + 37)] 1 ^ 37 . The requirement 7 < —2/3 asserts that X < 1. 

(v) Finally, when = we find x = 1 regardless of oj, 7 and /3, as expected. 

In Table 3 we illustrate the solutions to (|M]). They should be compared with the CDM with u> = 7, which are 
given by the first column in Table 2. They drastically differ for uj < —2/3. 



16 If x > 1, then 1^ 3> 1 for (3 > 1. Consequently, in this limit x = [— (1 + 3CJ)] 1 / 3 ". For consistency, the requirement x > 1 demands 
> -2/3. If S < 1, then i' 3 < 1 for > 1. Thus, in this limit 2 = [- (1 + 37)] 1/37 . The requirement x < 1 demands 7 < -2/3. 
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Table 3. x for various co, 7 and f3 3> 1 


CO 


7 = -0.4 


7 = -2/3 


7 = -5/6 


7 = -l 


7= -1.5 


-0.4 


3.824 


3.824 


0.850, 1", 3.824 


0.794, 1+, 3.824 


0.757, 1+, 3.824 


-0.5 


1.587 


1.587 


0.850, 1, 1.587 


0.794, 1+, 1.587 


0.757, 1+ 1.587 


-2/3 


1+ 


1 


0.850 


0.794 


0.757 


-1 


1" 


1" 


0.850 


0.794 


0.757 


-1.5 


1- 


1- 


0.850 


0.794 


0.757 



In terms of the cosmological redshift the solutions to (|90[) are 



[-(l + 3w)] 



l/3w 



n 



(AT) 



(DE) 



-1, 



l/3u 



(l + 3 7 )] 



l/3 7 



(92) 



Solution (i) generates models that coincide with a subset of the 7 = 10 models given by (|72|) where q crosses from 
plus to minus at z = z u . Solutions (ii)-(iv) show that an evolving DE with uj ^ 7 may affect not only the epoch at 
which q vanishes, but also the number of times it changes its sign. 

Models with co and 7 in the range given by (ii) have dW/dx > and q goes through zero three times: the first 
change is from plus to minus at z = z~ ( . The next change is from minus to plus at z — z*. Finally, there is another 

transition from plus to minus at z = z u . As an illustration, let us take co = —0.5, 7 = —2/3 and Qq M ^ — 0.3. From 
Table 1 we get xq = 1.759. Using the numbers given in Table 3 we get z 7 = 1.069, z* = 0.759, z u — 0.10. These 
solutions are similar to the one illustrated in Fig. [2] 

Models constructed from solutions (iii) have dW/dx < and q goes through zero only once, from plus to minus. 
The remarkable feature here is that the transition is at x 1 for all co < 2/3 and 7 > 2/3. However the observable 
redshift of transition z — z* does depend on co. 

Models generated by solution (iv) are interesting because the transition from positive to negative q depends only 



on 7, so it occurs at the same x for all to < —2/3. But, the corresponding redshift z 
co. 



does explicitly depend on 



6.2 Solution for 7 = 0: DE evolving from non-relativistic dust 



In (|9Tj) we have assumed 7^0. However models with 7 = are interesting because - following our discussion at the 
end of section (3.2) - they behave like pressureless dust at early times when x is small and as DE with 

p (DE) = wp (DB) 



(93) 



at late times. It turns out that (|89|) for 7 = admits a nice solution, viz., 

JDE) \ -P/ 3u 



x = 



(M) 



- 1 



Since the quantity inside the brackets must be positive we should require (See Figure [5]) 

nT < ' 



1 + ' 

At the current epoch the EoS and deceleration parameter are given by 



(94) 
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Figure 5: The decreasing functions represent the behavior of f2( M ' for models with 7 = and various choices of uj. 
They decrease from fiW = flmil = [l + 2 3w /' 3 ] at early times to f2( M ) = as x — !• 00 at late times. Clearly, 

1/2 < Vt^ax < 1 for any uj < and /3 > 0. Keeping w fixed we find that Qmi x — > (l/2) + for j8 3> 1 and £lmaL —> l - 
for p — > 0. In the case under consideration, for /3 = we recover the usual CDM model with W — uj/ '2. The 
increasing functions give 



1 

1- - 



n 



(AT) 



-,3/3^ 



n 



Thus, W -> w and q -> (1 + 3w) /2 as ^ DB) -> 1. 



Wo 



Chaplygin gas: If in addition to 7 — we set w = — 1 then the energy density and pressure of the DE component 
satisfy the EoS (|50l) . When (3 changes in the range 3 < /3 < 6 we recover the so-called generalized Chaplygin gas. For 
(3 = 6, which corresponds to the original Chapligin gas, the solution to is x = 1.069. To obtain specific numbers 
for the redshifts we need to feed the equations with f2[ ) M ' ) . The model requires VL M ^ < 2 — \/2 w 0.586. Thus, as an 
illustration we take tt ( M) = 0.3. With this choice we find xq = 1.465, Wq = —0.908 and qo = —0.453. The universe 
becomes dominated by the gas at z* = 0.465 and the phase of accelerated expansion starts at z = 0.370. 



6.3 Solution for uj = 0: dark energy transforming into non-relativistic cosmic dust 

Thus far we have assumed uj ^ 0. Models with u> — are diametrically opposite to the ones with 7 = 0. They 
behave as DE with p( DE > = jp( DE ) at early times and like pressureless dust at late times. Immediately after (|44[) . 
we mentioned that the replacement f3 — > —(3 changes the asymptotic role of uj and 7. The solution to ([89]) is now 
given by 
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Figure 6: The decreasing functions represent the behavior of 57^ M ' for models with w = and various choices of 7. 
The increasing functions give Q( DE \ As the universe expands fiW -> when uj ^ 0. However, il^ -> fi^jfj = 
[l + 2-^/^] -1 when cj = 0. Clearly, < Q^l < 1/2 for any 7 < and /3 > 0. Besides, ~> (V 2 )" for /3 » 1 

and fi^f^ — > for /3 — > 0. In the case under consideration, for j3 — we recover the usual CDM model with DE 
characterized by a constant EoS, viz., W = 7/2. 




Figure 7: The figure illustrates that accelerated expansion may occur in models with u = for different choices of 
7. The deceleration parameter has a minimum, say q m im around the transition epoch x = 1. For a fixed j3 the 
minimum is negative or positive depending on 7. In the Figure f3 — 10. For this value q m in is negative if 7 < —0.943. 
In general, the larger the magnitude of 7, the greater the magnitude of q m in ■ Also the accelerated expansion starts 
earlier and finishes later with the increase of |7|. 
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n 



(M) 



-/3/3 7 



-1//3 



X 



9. 



(DE) 



(95) 



which can readily be obtained from (|93j) by replacing j3 — > —(3 and uj — > 7. However, the analogy is not complete; 
namely instead of (|94[) - with the corresponding changes - in the case under consideration the condition on fig ^ is 
reversed, viz., 



1 



! + 2-37/^ ' 



Since /3 > 0, xq —¥ 00 as 57 



(Af) 



[l + 2 -37 /' 3 ] 1 and q — >• 1/2 as corresponds to dust models. In Figure [6] we 
illustrate this for various choices of 7. 

The question of whether or not this type of DE leads to an accelerated expansion depends on the choice of 7 and 
f3. This is illustrated in Figures [7] and |U 



Reversed Chaplygin gas: For the particular choice 7 = — 1 we have 



_ 



with A 



3 + /3 



and 



JDE) 



Co 



Co 



{it 



\(3+/8)/i8" 



For /3 > 0, the fluid behaves like a cosmological constant for x -C 1 and dust at late times. Matter dominates for 
small x so that q = 1/2 at early times and at late times. The expansion becomes accelerated in between for /3 > 9. 
See Figures El HI 



6.4 DE acting as a cosmological constant at late times: oj — — 1 

To obtain an explicit expression relating p^ DE ^ and p( DE ^ we use (|43|) to isolate x and then substitute the result into 
(|46|) . To ensure that the DE behaves as a cosmological constant at late times we set u = — 1. Solving for p( DE "> we 
obtain 



P 



(DE) 



P 



(DE) 



7 - (7+1) 



JDE) 



Co 



-0/3(7+1) 



7^-1. 



(96) 



This can be interpreted as the superposition of two non-interacting fluids, with energy densities p\ and P2, evolving 
according to the EoS 

Pi = IPX, p 2 = -V 1 - /3/3(7+1)] , (97) 
where A = (7+ 1) 7+1 ''. Thus, from the continuity equation we get 



Pi = 



2-3(1+7) : 



92 = 



A 



B 2 



xP/^+i) 



3(l+7)//3 



(98) 



where B\ and B2 are constants of integration. The above expression requires 7 > — 1, i.e., A > 0, otherwise pi 
diverges at some finite x ^ 0. For 7 > —1, p\ acts as quintessence and pi at early times - when x is small - behaves 
as pressureless dust: P2 ~ l/rr 3 . At late times /?2 behaves as a cosmological term since P2 = —P2 in that limit. For 
7 = and /3 = 6 we have a combination of dust and Chaplygin gas. 

It is well known that a scalar field (j> minimally coupled to gravity with a potential V(4>) can serve as a model 
for DE (see e.g. [T8]). In a DE-dominated universe a first order formalism gives {d4>/dxy = 3e (1 + W) /x 2 and 
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Figure 8: The figure illustrates the role of /?, which measures the rapidity of transition from u to 7. The larger 
j3 the shorter the transition. Comparison with Figure [7] illustrates that increasing j3 decreases the duration of the 
accelerated expansion and raises the magnitude of q m in- In this respect, the effects of 7 and (3 are very much alike. 
This is because, when u> = 0, both affect the slope of W in a similar way, viz., dW/dx ~ —7/3/4 near x = 1. 

V = (1 — W) p/2, where e = 1 or e = — 1 depending on whether <f> is an ordinary scalar field (W > —1) or a phantom 
field (W < -1). 

For lu = — 1 wc get 

1 



sinh 2 r]<f) ' 77 2V3VTT7' 



= ^ [(cosh ?7 0) 6(7+1)/ ' 3 + (1 - 7 sinh 2 7 7 0) (cosh7^) 2[3(7+1) ^ 1/ ' 3 
For 7 = and f3 = 6 this expression gives the potential for the Chaplygin gas discussed in [18] 



(99) 



6.5 DE acting as non-relativistic cosmological dust at late times: 7 = — 1 

Following the same steps leading to (|M|) . but this time setting 7 = — 1 we obtain 



p (DE) _ (DE) 



U - (Ul + 1) 



c 2 



(100) 



Once again the DE component can be interpreted as the superposition of two non-interacting fluids, viz., 

7 [l+0/3(w+l)] 



with A = (u + 1) C ; 



-/3/3(w+l) 



and 



Pi = 



gi 

a;3(n-w) ' 



P2 



So 



-1 -3(l+o>)//3 



2G 



where B\ and B2 are constants of integration, and A > to avoid singularities at some finite x. In this case the 
evolution of P2 is diametrically opposed to that discussed in Specifically, now P2 = —p2 = A~ 3, ( u+1 ^ ^ 

at early times and P2 ~ 1/x 3 at late times. A particular example of this model is the case where lj = 0, which 
corresponds to the reversed Chaplygin gas discussed in Section 5.3. 

Following the same procedure used in the preceding subsection, we can obtain a potential by treating DE with 
7 = —1 as an ordinary scalar field. We find 

^=sinh 2 # , V = ^=- 

and 

V{4>) = 9± [(cosh#r 6(1+aj)/ ' 3 + (1 - W sinh 2 #) (cosh#r 2[3(1+w)+,31//3 ] . (101) 

In order to avoid a possible misunderstanding, we should mention that the models with 7 = 0, discussed in 
Section 5.2, cannot be expressed as a superposition of fluids, except in the case where ui = —1. The same is true for 
the models with ui = of Section 5.3, which can be represented as a superposition of fluids only when 7 = —1. 



6.6 Parameter space 



In the CDM models there is a simple relationship between the cosmological parameters at the current epoch. This 
is given by the third equation in (IT2"j) . In the case under consideration, we can combine (j6"3")l . with n = 0, and 
to eliminate xq and obtain a single equation, viz., 



3( 7 -c)^ g) 
1 - 2q + 3 7 n ( DE) 



1 -7 



2qo — 1 — 3wf2 ' 
3 (7 - u) ti 



(DE) 



(DE) 



— 2( w -7) 



n 



(AO 



L fi 



(DE) 



/3/3 



(102) 



This replaces (|72|) when ^^7. 

Apart from this, some general constraints on the current values of the parameters can be obtained from ()63l) . 
Solving for Xq we get (w 7^ 7) 



3 7 ^ B) 



2?o + l 



3w 



(DE) 



-2q + l 



For £Iq MS> < 1/2 the condition xq > 1 gives: 



'0 

• For (oj - 7) > 0: 

Accelerated expansion requires 



1 1 3 / . \ n (DE) 



— w f2 



(DE) 



(iO + 7) < 



30. 



(103) 

(104) 
(105) 



Since cj > 7, it follows that 7 < — l/3fl 
• For (w - 7) < 0: 



(DE) 



1 3 ,, 

- + -UjQ 

2 2 



(D£) 



(106) 



The r.h.s. part of this inequality, together with Qq DE ^ > 1/2 and u < 7, leads to 7 > —2 (1 — 2qo) /3. 

• It should be noted that the inequalities (I104[) , (1106[) guarantee the positivity of the quantities inside the brackets 
in f|102[) for ui > 7 and w < 7, respectively. 
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7 Summary and conclusions 



Models in science perform two fundamentally different - but complementary - tasks. On the one hand, a model can 
indicate the type of phenomena that could actually occur in the context of the theory. On the other hand, it may 
shed light upon some problems of the theory. 

The CDM models with W = uj = constant are important first approximations in the description of DE. They 
satisfactorily explain the observational data for small values of z, which seems to be reasonable since the equations 
of state of the constituents of the universe have probably not changed much in recent times. The natural question to 
ask is whether these models will still be useful when we have access to observations of higher values of z. One would 
expect that the details of the expansion of the universe, which are determined by the evolution of its components, 
will become increasingly important as we look into the past. 

In this paper we have constructed a simple cosmological model which is physically reasonable, mathematically 
tractable, and extends the study of CDM models to the case where both the EoS for matter and DE vary with time. 
It is given by equations (IM)) - (I571) . The deviation from CDM models is measured by the functions f(z) and g(z) 
given by (|65[) . For CDM models f(z) — g(z) = 1. In our model f(z) pa 1, g(z) ps 1 for small z and f(z) pa (1 + z) 3n , 
g{z) ps (1 + 2 ;)~ 3( - a; ~ 7 ' ) for large values of z. Consequently, the model constructed here approximates the CDM ones 
for small z but significantly differ from them for large z. 

We have done a thorough review of all possible cases in the whole range of parameters of the EoS |43|) . (|53j) 
and discussed the corresponding physical models. Some of them are quite exotic as, e.g., those with uj > —2/3 and 
7 < —2/3 (see Table 3) that go from decelerated to accelerated expansion, and vice versa, several times. An example 
is illustrated in Fig [2j 

In section 5 we investigated the effects of a primordial EoS on the present accelerated expansion. In the CDM 
model there is a simple relationship between the current values of the cosmological parameters q, ui, QS DE > , which is 
given by (fT2")l . When ic/O that simple equation is replaced by three inequalities, viz., (|5T1) . (|5SJ), (|57|) for fi Q M ' < 1/2, 
fio = 1/2 and fijj > 1/2, respectively. They explicitly depend on n in such a way that setting n — we recover 
the CDM expression flTg) . 

We also studied the asymptotic model a 3> n which corresponds to the case where w abruptly changes from 
w = n to w — 0. For uj > —2/3 and Oq < 1/2 this and the CDM model are indistinguishable in the sense that 
they give the same values for xq, x, z* and z, regardless of n. The situation is quite different when uj < —2/3. 

We note that the solutions for x sharply depend on whether DE is phantom (uj < —1) or non-phantom (uj > — 1). 
Specifically, when the universe evolves from a radiation dominated early stage (n = 1/3) from the solutions (ii)-(v) 
of section 5.2 we find that x = 1 for — 1 < uj < —2/3 and x = [— (1 + 3w) /2] 1 ^ 3 " ^ for uj < —1. A similar situation 
is seen in where the limits on go are clearly dependent on the nature of DE. 

To illustrate the early and late behavior of the model, we consider the density parameter 

Q(M ) = (K + iy 3n/a (Kx* + lf n,a 

(K + iy 3n / a (Kx a + lf n/a + x 3(n- 7 ) 2 3(^- 7 )//3 (y + ^-BCw-rV/J ' 

Near the origin, for small values of x, there are three different cases: (i) fl^ 1 ^ = 1, (ii) fi( M ) < 1 and (hi) f}( M ) = 
corresponding to (n — 7) > 0, (n — 7) = and (71 — 7) < 0, respectively. The models discussed in section 5 have 
n > and u) = 7 < so they illustrate case (i). The models studied in section 6.2, which include the generalized 
Chaplygin gas, have n = 7 = so they illustrate case (ii). We have restricted our discussion to n > 0, 7 < so 
case (iii) has been excluded from consideration. In the limit x — > 00 there are two possibilities: Q( M *> — > and 
1/2 > fj( M ) > for uj < and ui = 0, respectively (we have not considered ui > 0). The models studied in section 
6.3 illustrate the second case. 

In this work we have not discussed the physical mechanism leading to DE models that cross the W = — 1 barrier. 
This is an important and interesting question but this is out of the scope of the present work. 
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